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We investigate the simplest quintessence dissipative dark matter attractor cosmology character- 
ized by a constant quintessence baryotropic index and a power-law expansion. We show a class of 
accelerated coincidence-solving attractor solutions converging to this asymptotic behavior. Despite 
its simplicity, such a "superattractor" regime provides a model of the recent universe that also ex- 
hibits an excellent fit to supernovae luminosity observations and no age conflict. Our best fit gives 
a — 1.71 ± 0.29 for the power-law exponent. We calculate for this regime the evolution of density 
and entropy perturbations. 
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I. INTRODUCTION 



A new picture of the universe is emerging from observations of large scale structure, searches for Type la supernovae, 
and measurements of the cosmic microwave background anisotropy. They all suggest that the universe is undergoing 
cosmic acceleration and is dominated by a smoothly distributed dark energy component with large negative pressure 
jl] Ej|. Moreover, recent data coming from BOOMERANG and MAXIMA projects seem to indicate that the density 
of the universe is near the critical density The most frequently proposed candidates for this dark energy are a 
cosmological constant, or vacuum energy density [^), and quintessence Q, a scalar field with negative pressure. On the 
other hand, gravitationally clustered matter is usually assumed to be dominated by cold, collisionless, nonbaryonic 
dark matter (CDM). 

In addition to the old cosmological constant problem, related to the smallness of the observational upper bound on 
the vacuum energy density compared to particle physics scales jjj, a new challenge to the model with cosmological 
constant and CDM (ACDM) is the "cosmic coincidence" problem: why is it that the vacuum density dominates 
the universe only recently S. A dynamical self-interacting scalar field may explain why the dark energy is small. 
However, a minimally coupled scalar field combined with perfect fluid matter, such as in "tracking" QCDM models 
cannot explain the observed acceleration and solve the coincidence problem pO|Jll] ]. On the other hand it has 
been shown in Ref. ]l0[ | that both acceleration and coincidence can be satisfactorily explained by a combination of 
quintessence and dissipative dark matter (QDDM). For these models it was shown that late-time attractor solutions 
exist with very interesting properties: an accelerated expansion, spatially flatness, and a fixed ratio of quintessence to 
dark matter energy density. Recently, it has been shown that such accelerated coincidence-solving attractor solutions 
also exist in spatially flat models with an exponential self-interaction quintessence potential and phenomenologically 
chosen couplings between the quintessence field and baryotropic dark matter Jl^,|l3[ . 

Consideration of dissipative effects in dark matter also arises from increasing evidence that numerical simulations of 
dark matter halos on sub-galactic scales based on conventional CDM models lead to conflicts with observations. One of 
the major problems shown by these simulations is that galactic halos are too centrally concentrated fll4[ . Confirmation 
of this problem would imply that structure formation is somehow suppressed on small scales. Several scenarios 
addressing this issue have been considered assuming some kind of interaction for dark matter particles: nonthermally 
produced and weakly interacting jL5|, self-interacting [^6|, repulsive ]l7| , annihilating |l8[| , and decaying [jl9| . It is 
quite reasonable to expect that dark matter is out of thermodynamical equilibrium and these same interactions are 
at the origin of the cosmological dissipative pressure. A simple estimation shows that a cross section of the order of 
magnitude proposed in these halo formation scenarios, corresponding to a mean free path in the range 1 kpc to 1 
Mpc, yields at cosmological densities a mean free path a bit lower than the Hubble distance. Hence a description for 
interacting dark matter as a dissipative fluid is valid at cosmological scales |2(i[| . 

Dark matter may also have several components, from very heavy weakly interacting massive particles to a lightweight 
neutrino El). Even if this neutrino- like component does not contribute significantly to the density budget of the 
universe 3|7 it may have nevertheless a relevant dynamical role as distinct components of dark matter, cooling at 
different rates, give rise to bulk viscosity p2| . 

Recently the relationship between dark matter clustering, the nature of dark matter and the origin of ultrahigh 
energy cosmic rays (UHECRs) has been explored §§. Heavy particles (m x ~ 10 12 - 10 14 GeV) can be produced in 
the early universe in different ways [p4|-p7| and their lifetime can be finite though longer than the age of the universe. 
In these circumstances, super-heavy particles can represent an appreciable fraction of dark matter and the decay of 
these particles results in the production of UHECRs, as widely discussed in the literature |2^ , |28| -^o|] . In particular, 
if the relics cluster in galactic halos, as is expected, they can explain the cosmic ray observations above ~ 5 x 10 19 
eV. The effect of decaying dark matter on the equation of state has been studied numerically in Ref. J3l| , and the 
equivalence between particle production and dissipative bulk viscosity has been investigated in Ref. Q . 

All this shows that many different scenarios may occur where significant dissipative processes develop in dark 
matter, in particular when it behaves as a viscous fluid. So a variety of accelerated coincidence-solving attractors 
are possible. In this paper we will show that for a wide class of QDDM models attractor solutions are themselves 
attracted towards a common asymptotic behavior, the "superattractor" . The superattractor scenario is described 
in Sec. II and observational constraints on this scenario are investigated in Sec. III. These include the luminosity 
distance-redshift relation for type la supernovae and the age of the universe. We apply the covariant gauge-invariant 
formalism to calculate the evolution of density and entropy scalar perturbations in the QDDM regime and solve 
them for the superattractor regime in Sec. IV. Finally the main conclusions are discussed in Sec. V. Units in which 
c = 8irG = ks = 1 are used throughout. 
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II. QDDM SUPERATTRACTOR SCENARIO 



In Ref. it was shown that the Friedmann-Lamaitre-Robertson-Walker (FLRW) universe filled with perfect 
normal matter plus quintessence fluid, corresponding to some minimally coupled scalar field governed by the Klein 
Gordon equation, cannot at the same time drive an accelerated expansion and solve the coincidence problem. To 
solve it, some additional bulk dissipative pressure ir in the stress-energy tensor of dark matter was considered. Any 
dissipation in FLRW universes has to be scalar in nature, and in principle it may be modeled as a bulk viscosity effect 
within a nonequilibrium thermodynamic theory such as the Israel-Stewart theory fl33| , [34| . In a certain regime, that 
formulation can be approximated by the more manageable truncated transport equation 

7r+T7T= -3C-H", (I) 

where H = a/a denotes the Hubble factor, ( stands for the phcnomenological coefficient of bulk viscosity, and r is the 
relaxation time associated with the dissipative pressure [p5|-p7[. As usual an overdot means derivative with respect 
to cosmic time. 

The overall stress-energy tensor of the QDDM model reads 

Tij = {p m + P4>+ p m + + Tr)uiUj + {p m + p0 + ir)gij , (2) 

where p = p m + p^, and p = p m + p^. Here p m and p m are the energy density and pressure of the matter whose 
equation of state is p m — ("j m — l)p m with baryotropic index in the interval 1 < j m < 2. Likewise p^ and p<p, the 
energy density and pressure of the minimally coupled self-interacting quintessence field <p, are related by the equation 
of state, p$ — (70 — l)p0, with baryotropic index 

where for non-negative potentials V(<f>) one has < 7^ < 2. The scalar field can be properly interpreted as quintessence 
provided 70 < 1 -see e.g. Q. In general 7^ varies as the universe expands, and the same is true for j m since the 
massive and massless components of the matter fluid redshift at different rates. 

The Friedmann equation and the energy conservation of the normal matter fluid and quintessence (Klein-Gordon 
equation) are 



k 1 
3 ( 



H 2 + — = -(p m +Pcb) (k = 1,0,-1), (4) 



p m + 3H [ 7 m + — ) p m = 0, (5) 
V Pm ) 

P4, + 3ff 70/90 = 0, (6) 

where the prime means derivative with respect to <fi. Introducing f2 m = p m /Pc f^,= Pc/>/pc, with p c = 3H 2 the 

critical density and flk = —k/(aH) 2 plus the definition = ft m + Q^, the set of equations can be recast as 
(cf. i|) 

n m + n 4> + n k = 1, (7) 

Q m = , (8) 




tl^ + ZH U^ + ^L ) fi = O, (9) 



3H 2 



where 7 is the effective baryotropic index given by 
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7^ = 7m^m + 70^0- (10) 

Equations ©-© have fixed point solutions Q = 1, Q m = fl ma and £1$ — fl^a, respectively, when the partial 
baryotropic indices and the dissipative pressure are related by 

7T 2H 

7m + — =70 = -T7T72- (11) 
Pm °- rz 

Asymptotical stability of f2 = 1 occurs whenever 7 + n/ p < 2/3. This condition, together with Eq. (fnl), leads to the 
additional constraints it < (2/3 — 7 m ) p m , which is negative for ordinary matter fluids, and 70 < 2/3. Additionally 
the stability of Q m — fl ma and f2^ = 0^ a , and hence of the ratio i} m /fl ( j ) , has been established in p0| . These solutions 
provide a natural explanation to several features observed in our universe: an accelerated expansion, spatial flatness, 
and a ratio of dark energy to matter density of order unity. We denote by a subindex a the asymptotic limit of 
magnitudes in the attractor regime, while the subindex will denote as usual their current values. 

Combining Eq. ([ll]) with Eq. (jl|) we obtain the equation of motion of the attractor solutions of the system (0) , 
@> ©j © satisfying flatness, acceleration and coincidence: 

,-(1 + 3^)^ + ^-^ = 0. (12) 

Here v = (tH) 1 is the number of relaxation times in a Hubble time - for a quasistatic expansion v is proportional 
to the number of particle interactions in a Hubble time. Perfect fluid behavior occurs in the limit v — > 00, and a 
consistent hydrodynamical description of the fluids requires v > 1. Rewriting (|l^) in terms of the field baryotropic 
index 70, we get 

70 = 370 - + 37m)70 + v (j m - jj ) > (13) 

where a prime indicates derivative with respect to rj = In a. 

When the phenomenological coefficient of bulk viscosity satisfies ( = £ s , where 

C, = Q m a (7m - 70s) [l - 370s v~ X ] H = kH , (14) 

Eq. ( |l3| ) admits the constant solution 70 = 70 S . It gives an accelerated expansion in the late time regime when 
70s < 2/3. As C > and 7 m > 1, the hydrodynamical parameter v is restricted to v > 370 S . The case of constant k 
in the interval < k < 1 arises, for instance, in a radiating fluid, and the nearly linear regime, with slowly varying v 
and 7 m , was already investigated in the quasiperfect limit, corresponding to v^ 1 — * J1Q|. 



To analyze the stability of the solution 70 = 70 S we insert Eq. (|14j) into (13) 



70 = 3 (70 - 70s) ~{v + 37m) (70 - 70s) • (15) 



As 70 S < 2/3, v > max(370 S , 1), and j m > 1, Eq. ( J15[ ) shows that dj'^/d"/^ < in a neighborhood of 70 S . Hence 
this constant solution is asymptotically stable, showing that all solutions of Eq. (|l2|), that is, all the accelerated 
coincidence-solving attractors of the system (|]), (||), (||), (|l|), are themselves attracted towards the constant solution 
70 — 70s provided that they satisfy £ ~ Cs when t — * 00. As the same occurs with all solutions whose initial conditions 
lay within the domain of attraction of each of these attractor solutions, we will refer to the constant solution as the 
"superattractor" of this class of QDDM models. We denote with subindex s magnitudes in the superattractor regime. 
As an example of models satisfying £ — > ( s we note the case £ cx ^/p m , investigated in Refs. (3^-Q. 

Let us characterize this asymptotic stage. From Eq. ([ll]) the superattractor solution yields a power-law evolution 
for the scale factor 

a s (t) = a (j-^J , (16) 

where a — 2/370 S , ao is the current scale factor and t s is the age of a superattractor universe. The dynamics of the 
scalar field in the attractor regime is obtained from Eqs. (||) and (pT|): 

V = ^ a (2- 7(t> )H 2 , (17) 
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fa = 3^70 H . (18) 

These equations together with (12) close the problem of finding the potential V and the scalar field <f> as functions 
of t. Expressing the potential (17) and field (fijl) in terms of cosmological time, that is, V = f20 a (3-ff 2 + H) and 



,2 



-2n,p a H, it is easy to find, for the superattractor solution, 



W) = r^3(2-7*.)e- A *, (19) 
3 {J4,sts) 



*-»-*> + V^4' ,20) 

where A — (370 S /fi0 a ) 1 / 2 is the slope parameter. Then we find a = 2 / (Q^aA 2 ) and fa = (2 /A) hi70 S . We note that 
in the superattractor regime, the exponent a is larger, by a factor of l/Q^a, than the exponent 2/ A 2 of a scalar field 
dominated era. This is a particular instance of a general property of QDDM models in that dissipative effects assist 



quintessence driven acceleration through the negative non-equilibrium pressure tt. Using Eq. (11) we can evaluate its 
ratio to quintessence pressure: 

_ 70 - 7m ^ma 

Hence, in the superattractor regime, this ratio is also a constant (provided that j m is a constant), and when dark 
matter is cold it is just the ratio of matter to quintessence energy density. 

Let us sketch the evolution of the actual universe from a nearly thermodynamical equilibrium early era (when 
ir/p m ~ 0) into this superattractor stage (when ir/p m — 70 S — 7 m ). Along the radiation and matter dominated eras 
\^l Pm\ <C 1 and the inequality j m + ir/p m > 70 s holds. Then dissipative processes become more significant, the 
attractor condition ( |lT| ) is approached, 70 is driven towards 70 s , and the expansion of the scale factor accelerates. 
Clearly the timescale of this transition period depends on the details of the dissipative processes occurring in dark 
matter, encoded into the evolution of the dissipative magnitudes and v(t) - some models exhibiting this con- 
vergence stage have been investigated in Ref. |Hj. The superattractor stage finally settles when C — Cs- Assuming 
that the ratio of this transition period to the age of the universe to is small and this transition period started early 
enough, we may approximate the recent evolution of a(t) by the superattractor solution (|l(]); hence, 700 — 70 S and 
(1 + z)- 1 = a/ao- {t/t ) a - 

As Eq. (|19| ) shows, the class of models converging to the superattractor stage has V(fa ~ V s (fa for large <j>. In 
addition to this exponential tail, no other constraint has to be imposed on the quintessence potential for convergence 
to the superattractor era. For a wide range of initial values of cj> and <fi the quintessence field approaches a common 
evolutionary path j20| ) for which 7^ = 70 S ; i.e., the late behavior is insensitive to the initial conditions. 



III. OBSERVATIONAL CONSTRAINTS 
A. Luminosity distance of supernovae 

It has been found that supernovae of type la (SNela) are nearly standard candles. Properly corrected, the difference 
in their apparent magnitudes is only related to differences in luminosity distance dh and consequently to cosmological 
parameters . Taking profit of this property, several accelerated expanding cosmological models like ACDM and 

QCDM have been fitted to recent observations of high redshift supernovae {z<\) |l3U45| Though they have a 
good fit in some regions of the parameter space corresponding to an accelerated expansion, these models require fine 
tuning to account for the observed ratio between dark energy and clustered matter. 

On the other hand, QDDM models provide models that simultaneously provide an accelerated expansion and 
solve the coincidence problem within the general trend of the universe towards an attractor and with it towards the 
superattractor. As in QCDM models, these scenarios depend on the quintessence potential and initial conditions. In 
addition, they depend on the evolution of the magnitudes characterizing the details of dissipative processes occurring 
in dark matter. Let us examine the issue of reconstructing part of this evolution through observations of distant 
SNela. 

Ignoring gravitational lensing effects, the standard expression for the luminosity distance to an object at redshift z 
in a spatially homogeneous and isotropic universe is p8 
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d L (z) 



l + z 



dz' 



ffol^fcol 1 / 2 " V Jo H(z'), 

with S(u) = (sin it, u, sinhu) for k = (1,0,-1) respectively. Then for a QDDM universe we have 



(22) 



H(z) 



f2 m0 exp 



dz' 



7m + 



1/2 



l + z' 

+Sl 40 expUj Jf—^+n k0 (l + z') 2 ^ (23) 
In the case of a spatially flat universe, dh has a simpler expression in terms of the effective baryotropic index 



dh{z) 



1 + z f z , , 3 [*' dz" 

/ dz cxp / 

Ho Jo 2 Jo 1 + 2" 



7 + 



where 



7T 

7+- 
P 



7T 

7m H 

Pro 



7m H 70 I ^00 exp I 3 

Pm J \ JO 



dz' 



l + z 



ll<t> 



(24) 



(25) 



We see from Eqs. (|24|) and ( p5| ) that di,(z) depends on the time evolution of dissipative processes through a double 
integral in n/p rn (z) and on the time evolution of the quintessence field through a double or triple integral in 7</,(z) (we 
assume that 7 TO ~ 1). On one the hand, there is a degenerancy here as these two functions cannot be reconstructed 
from knowledge of the single function d^jz). Besides the time variation of these magnitudes is largely smoothed 
out and, similarly to the analysis in Refs. |49|-|5l|| for QCDM models, we find that the luminosity distance is highly 
insensitive to these variations. So we may safely replace these time-varying magnitudes by their mean values in the 
interval (0, z). Assuming that the universe has already settled in the superattractor regime at the age of the farthest 
observed supernova the degenerancy is eliminated and both functions become constant. Then, the expression ( p? 
simplifies drastically and we obtain (cf. |32| ) 



dh{z) 



(l + z) [(l + z f-l] 
PH 



(26) 



where (3 = ad /a 2 = 1 — 1/a is the acceleration parameter. So (3 increases with a and 1 < a < 00 corresponds to 
< j3 < 1. On the superattractor we have 



70s 



2 

3a 



:a-/3). 



(27) 



We have used the sample of 38 high redshift (0.18 < z < 0.83) supernovae of Ref. p3j, supplemented with 16 low 
redshift (z < 0.1) supernovae from the Calan/Tololo Supernova Survey |53]| . This is described as the "primary fit" 
or fit C in Ref. where, for each supernova, its redshift Zj, the corrected magnitude rrii and its dispersion <7i were 
computed. 

The predicted magnitude for an object at redshift z with luminosity distance d^ is 



m(z) = M + 5\ogV L (z), 
where M. is related to the absolute magnitude M by 

'ifopcm/s Mpc] 



M = M - 5 log 



;[km/s] 
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(28) 



(29) 



and "Dl is the luminosity distance in units of the Hubble radius. 

We have determined the optimum fit of the superattractor model by minimizing a x 2 function: 

2 _ i m i - m(z t ; P,M)} 2 



(30) 
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where N — 54 for this data set. The most likely parameters are (/3, M) = (0.395,23.96), yielding Xmin/^DF = 1-12 
(Ndf = 52), and a goodness-of-fit P(x 2 > Xmin) = 0.253. These numbers show that the fit of superattractor QDDM 
cosmology to this data set is as good as the fit of the ACDM model (see Fig. 1). We note that it occurs even though 
the superattractor model has basically only the acceleration parameter to fit large rcdshift supernovae (M. being 
largely determined by low redshift supernovae) . It also means that the density of clustered matter is not constrained 
by measurements of SNela and it has to be determined through independent observations. 

We estimate the probability density distribution of the parameters by evaluation of the normalized likelihood |54|] 

Then we obtain the probability density distribution for f3 marginalizing p(j3,M) over M. This probability density 
distribution p((3) is plotted in Fig. 2 and it yields (3 = 0.398 ± 0.104 (la). Hence we can state that 0.085 < < 0.711 
with a confidence level of 0.997, so that an accelerated superattractor QDDM universe is strongly supported by this 
data set, in agreement with a similar analysis of ACDM and QCDM models ||43| , [i~4] . p : 5|j46| . 



B. Age of the Universe 



For a QDDM model the age of the universe has the integral representation 

ta = Jo (! + *)#(*)' (32) 

where H(z) is given by Eq. (|2"5|). Along the transition period of the universe from an early nearly equilibrium 
stage towards the superattractor stage the inequality j m + ir/p m > holds. As a consequence H(z) > H s (z) — 
i?o(l + z) 1 ~' 3 (whenever fifc > 0) and to < t s = a/Ho. The characteristics of this transition period, and hence the 
value of the difference t s — to, are model dependent and will not be dealt with here. 

Using the fit to SNe la of the previous subsection and Eq. ( p7|) we obtain a = 1.711 ± 0.288. Combined with 
Hq = 65 ± 5 km/s Mpc ||, it yields for the age of the superattractor universe t s = 25.9 ± 4.8 Gyr. This shows that 
the QDDM cosmology does not suffer of any age discrepancy and can accommodate comfortably the la interval 9-16 
Gyr for the age estimate of globular clusters ||] . 



C. Parameters of the superattractor era 

Using the probability distribution for /3 and Eq. (|27]) we obtain 7^ s = 0.401 ± 0.069, in agreement with previous 
results for QCDM models in the limit f2 m — > 0. Then, using Eq. (|lTJ) and assuming that dark matter is cold 
(7m = 1), we find n/p m — —0.599 ± 0.069. This figure implies that substantial dissipative processes are taking place 
in dark matter. This fact is also shown by the large value of the effective baryotropic index. In effect, combining the 
estimate fl m o = 0.35 ± 0.07 from cluster baryons Q with the a priori constraint fl = 1 and inserting into Eq. ( |l0| ) 
we get 7 S = 0.61 ± 0.09. For a perfect fluid this value would correspond to a power-law exponent 2/37 s ~ 1.1, quite 
lower than a. Using Eq. (p"4]), we find that the linear relationship n = n\ + k^v^ 1 holds, where n\ = 0.21 ± 0.05 and 
k 2 = — 0.25±0.07. The requirement of asymptotic stability of the superattractor solution imposes that v > 1.20±0.21. 

Using Eqs. (|l7|) and ([18|), we find the quintessence kinetic energy density parameter flxo — 0.13 ± 0.03 and 
flvo — 0.52 ±0.06 for the potential energy density parameter. These figures show that the scalar field is moving down 
the potential outside the slow-roll regime. Similarly we find for the slope of the exponential potential A — 1.36 ± 0.14 
and for the current value of the scalar field cf>o = —1.3 ± 0.3. This implies a mass parameter of the Planck scale. 



IV. DENSITY PERTURBATIONS 



Besides the model degenerancy of luminosity distance determinations, even within the superattractor cosmology 
there are some thermodynamical parameters like v and the speed of sound c s that are not fixed. For this reason we will 
investigate the evolution of density fluctuations in the perturbative long-wavelength regime during the sup erattractor 
era. It is possible that weak lensing techniques could yield more information about these parameters ]5q |. 

Scalar perturbations are covariantly and gauge-invariantly characterized by the spatial gradients of scalars. Density 
inhomogeneities are described by the comoving fractional density gradient |57[] 
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Si = ^ , (33) 
P 

where stands for the covariant spatial derivative DjAi... = hj k hi l ■ ■ ■ VfeA;.... The scalar part 5 = aD l 8i = (aD) 2 p/p 
corresponds to the usual gauge-invariant density perturbation scalar e m |38| , |59| , which encodes the total scalar con- 
tribution to density inhomogeneities. Also the comoving expansion gradient, the normalized pressure gradient, and 
normalized entropy gradient are defined by p7| , p0| 

t = aDiti , pi = , e.i = , (34) 

P P 

n being the particle number density, T the temperature, and s the specific entropy per particle. The evolution equation 
for scalar density perturbations reads J(3(| 

6 + H (8 - 67 + 3c s 2 ) 6 - \R 2 {l + 5 (7 - i f - 6c s 2 + 
(1 - 3 (7 - l) 2 + 2c 2 ] kX S - c 2 s B 2 S = S[e] + S[tt] + S[q] + S[a] , (35) 



where 

' dp\ 1 /dp 



c 2 = I — I • r = 



dp J s ' nT \ds 



(36) 



are, respectively, the adiabatic speed of sound and a non-baryotropic index. The sources in the right-hand side of Eq. 
( |35| ) arising, respectively, from entropy perturbations, bulk viscous stress, energy flux, and shear viscous stress are 
given in |60[. Since in our case there are no shear viscous stress (cry = 0) and S[q] vanishes by choosing the energy 
frame (c/j = 0), we reproduce here only the expressions for S[e] and S[tt]: 

S[e] = r (3fcff 2 + D 2 ) e , (37) 
S[vr] = - (3kH 2 +T) 2 )B, (38) 



where the scalar entropy perturbation 



and the dimensionlcss perturbation scalar 



= aD J ei = ^^D 2 s (39) 



(40. 

p 

related to the inhomogeneous part of the bulk viscous stress, were defined. 
Also, the entropy perturbation equation in the energy frame is 

e + 3H (c 2 - 7 + 1 + r) e = -3HB . (41) 

The coupled system that governs scalar dissipative perturbations in the general case is given by the density pertur- 
bation equation (|35|), the entropy perturbation equation (|4l|), the equation for the scalar bulk viscosity (|IJ), and the 
equation for temperature perturbations. 

When only bulk viscous stress dissipation is present, the coupled system can be reduced to a pair of coupled 
equations in S (third order in time) and e (second order in time). For a flat background, the equations are |6C|]: 



r 5 + [I + 3 (2 - 7 + c 2 ) tH] 6 + H {8 - 67 + 3c 2 + 3r (c; 
-\ [-14 + 75 7 - 48 7 2 + (21 7 - 30)c s 2 ] tH } 6 
-\H 2 {6 - 10 7 + 5 7 2 - 6c s 2 - 4r (c s 2 )' 
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-2 [-6 + I87 ~ !57 2 + 57 3 + (6 - 28 7 + IO7 2 ) c 2 ] tH} 6 



fl2 Cn2 ^n2A , 3a2 (7 ~ 1)^2 f~2 X S , 2^2 



^D 2 (d 2 *) 



P7 



D 2 (D 2 6) + tc 2 D 2 5 + rrD 2 e 



(l-37rF)c 2 + r(c 2 )-+3(|^ 



D 2 <5 



and 



+ 



'l-37ri/)r + rr + 3 



9( 



D 2 e. 



re + [1 - | (-2 + 3 7 - 2c 2 - 2r) riJ] e 



(42) 



-3# 7 - 1 - c s 2 - r + 37 (7 - c s 2 ) ri? + r (c 2 + r) 



P fdC 



7 \<9s 



7 



7 



5. 



(43) 



We shall consider here the evolution of the density and entropy perturbations in the superattractor stage with the 
conditions r = and dv/ds — 0. Together with Eq. (|l4|) they imply 



/ 5C 



«1 + K 2 /v 



(44) 
(45) 



In this case, Eq. (f42|) decouples to give 



Cl 



C 2 c , c 3 



5 +—5 + ^5 + ^6 = c 4 i 2Q D 4 <5 + c 5 t 2a D 4 5 + c 2 s V z 8 



±2a r>4i 



f 2 



t 3 



c 7 D 2 <5 



and Eq. (43) becomes 



.. c s . Cg C10 1 Cn 



(46) 



(47) 



where the constant coefficients ci . . . cn depend upon the parameters of the model: v, Ki, K2, TO , a, 7, c 2 , and the 
present value of the scale factor ao- For our purposes, only the explicit expression for ci, cs, and C9 are relevant, being 

ci = av + 3a (2 — 7 — c 2 ) , 



c 8 = a 



^(37"2-2c 2 ) 



c 9 = -3a 2 [v (7 - 1 - c s 2 ) + 3 7 (7 - c s 2 )] 



(48) 



We deal with the system (|46|),(47) by performing separation of variables in the form 5 = S x S t and e = e x e t , where 
8 X and e x depend upon the spatial variables while St and et are functions of the coordinate time t. Then, Eq. ( |4^ ) 
can be recasted as 








(49) 

\ ut t I o x 

which can only hold if 

(D 2 - M ) 4 = , (50) 
/i being an arbitrary constant. Also, Eq. ( p7| ) leads to 

i 2 (e t /e t ) + c 8 * (e t /e t ) + c 9 4 



ClO 



Ut/etJ + cu (S t /e t ) 



(51) 



which requires e x = A 8 X , with A an arbitrary constant which can be absorbed into the temporal functions, resulting 
in the same spatial distribution of the entropy and density perturbations. Under the condition (|5(]) the evolution 
equation for St becomes 



Cl X i (° 2 2 ,2a 2\ r 



.f 3 



- M 2 i 2 « c5-/i(y + C7 )" 



<5 t = 



(52) 



As we are interested in the asymptotic behavior for the superattractor regime, it suffices to consider the dominant 
terms in (p2), a being a positive number, we have 



Cl 



St +j5 t -v 2 t 2a c 4 S t - ix 1 t 2a c 5 S t 







(53) 



Equations (Em and ( p3[) give the form of the asymptotic density perturbations in the model. The parameter fi 
appearing in ( |50| ) depends, in principle, on the boundary conditions of the problem being the quantity /i -1 / 2 , a 
characteristic coordinate length related to the range of the exponentially decaying modes of the spatial part S x . In 
the special case y, = 0, Eq. ( |50| ) has the Laplace form, describing long-range density perturbation modes. Here, we 
are going to study the asymptotic evolution of the long-range modes by performing a series expansion of St in powers 



of fi 2 . Up to first order we have St = S^ + fi 2 S^ 1 ' , then, replacing this expression in Eq. (53) and retaining terms up 
to first order in [i 2 we obtain 



,2 XW 



t 2a c 5 S 



(0) 



,(0) Cl " (0) 



The zeroth-order solution gives 



i2-ci 



(l-ci)(2-ci) 



t 



A 2 t + A 3 , 



sr = o . 



for ci ^ {1, 2} and Ai , i = 1, 2, 3 integration constants. Then, 5^ satisfies the inhomogeneous equation 



"° C 4 Si 0) 



t 2a c 5 <5 t (0) = 



(54) 



(55) 



(56) 



whose general solution has the form 



* t (1) = 



(1- C1 )(2- C1 ) 

+a 2 t 2a + 3 + a 3 t 2a+5 ^ 



B 2 t + B 3 + ai t 



2ct+4-ci 



04 t 



2a+4 



(57) 



where the coefficients Bi, i — 1,2,3 are integration constants depending upon the initial conditions, and ai, i — 1, . . . , 4 
depend on a, ci, and the Bi. Then, the asymptotic evolution of the long range modes up to first order in \i 2 results in 
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a combination of powers of time. When c\ < 2 all the exponents are positive, while for c\ > 2 there exists a decaying 
mode. With the parameters of the model, because of Eq. ( [48|) this last situation is equivalent to c 2 < v/Z + 0.999 
which is always satisfied because of the lower bound v > 1.2 (in our units c 2 < 1). As a consequence, a decaying 



mode exists in this model and the dominant exponent is 2a 

P(k,z) (x (l + z)-( 4+8 / Q ). 



On the other hand, Eq. ( |43[) becomes, in the leading regime, 



4. Hence the density power spectrum redshifts as 



Cs . Cg 2 
e-t + —e t + -j e t = exit 



1+- 7 -l 

v 



(58) 



whose solution is 



e*(0 



Cl3* 



2a+3 



(59) 



where Ai,2 are the roots of the equation A 2 + (eg — 1)A + eg = 0, B^^B^ are arbitrary integration constants, and 



Ci2,ci3 are functions of the parameters and the previously defined integration constants, 
perturbations may yield additional information about the parameters v, c 2 , and rema- 



in this way the entropy 



V. DISCUSSION 



We have shown that within the class of accelerated coincidence-solving attractor solutions of QDDM models, a 
distinguished attractor solution exists with constant quintessence baryotropic index. This superattractor cosmology 
is quite attractive because of its simplicity: the scale factor expansion follows a power-law, and all its parameters can 
be evaluated in a model independent way. Notwithstanding its simplicity, its fit to SNela observations is as good as 
the ACDM model, and it does not suffer from any age conflict either. 

Our results are quite valuable to give general statements about a large class of QDDM models that are driven 
towards this superattractor but they are limited to the late-time regime. Undoubtedly a large variety of models 
arise both from different quintessence potentials as well as from different dissipative processes in dark matter. The 
only requirement is that the potential has the exponential tail (|l9| ) or, equivalently, that the viscosity coefficient 
has the asymptotic behavior (p"4|). This diversity of possibilities makes the transition period of the universe from 
its nearly thermodynamical equilibrium early stage towards the superattractor regime model dependent and requires 
more detailed investigation. In particular it would be of importance to evaluate the effects of dissipative processes on 
the CMB angular power spectrum. 

Standard CDM models have been quite successful in describing the evolution of cosmic structure. Hence a suffi- 
ciently long matter-dominated era must have taken place during which the observed structure grew from the density 
fluctuations measured by CMB anisotropy experiments. As a consequence the transition to the currently observed 
accelerated regime should have been quite recent in the evolution of the universe. For QDDM models this seems 
to imply that dissipative effects in dark matter were small until density inhomogencitics became large. If so, it is 
suggestive to think that the size of dissipative processes, as measured by the ratio n/p m , grows with dark matter 
density and became large precisely because of the development of inhomogeneities. The observed smoothness of halo 
structure and the correlation between high energy cosmic ray production and dark matter clustering might provide 
support for this relationship. 

Even though the supernovae search is extended to z > 1 this does not enable a precise determination of the time 
variation in the quintessence baryotropic index 7^ and the ratio it / p m because the luminosity distance depends on 
these magnitudes through a multiple-integral relation that smears out detailed information about their variability. 
Besides their independent variation cannot be disentangled from the determination of a single function. Hence 
further independent cosmological probes are required to investigate this issue. As a first step in this direction we have 
calculated the dominant large-time long-wavelength behavior of density and entropy fluctuations in the superattractor 
regime. 

The combined distribution of dark mass and dark energy can be investigated via its gravitational effects inducing 
correlated shear in the images of distant galaxies. This weak lensing effect could, in principle, be used to probe 
the large-scale structure and thereby yield additional information about the thermodynamical and cosmological 
parameters not constrained by the luminosity-redshift relationship. 



11 



ACKNOWLEDGMENTS 

This work has been supported by the University of Buenos Aires under project TX-93. 



NA. Bahcall et al, Science 284, 1481 (1999). 

M.S. Turner, Phys. Rep. 333-334, 619 (2000). 

J.R. Primack, Nucl. Phys. B (Proc. Suppl.) 87, 3 (2000). 



A.H. Jaffe et al., Phys. Rev. Lett, (to be published), istro-ph/0007333 



M.S. Turner, in Critical Dialogues in Cosmology, edited by N. Turok (World Scientific, Singapore, 1997). 
R.R. Caldwell, R. Dave and P.J. Steinhardt, Phys. Rev. Lett. 80, 1582 (1998). 
S. Weinberg, Rev. Mod. Phys. 61, 1 (1989). 

P. Steinhardt, in Critical Problems in Physics, edited by V.L. Fitch and D.R. Marlow (Princeton University Press, Prince- 
ton, 1997). 

I. Zlatev, L. Wang and P.J. Steinhardt, Phys. Rev. Lett. 82, 896 (1999); I. Zlatev and P.J. Steinhardt, Phys. Lett. B 459, 
570 (1999); J.P Steinhardt, L. Wang and I. Zlatev, Phys. Rev. D 59, 123504 (1999). 
L.P Chimento, A.S. Jakubi, and D. Pavon, Phys. Rev. D 62, 063508 (2000). 

S. Weinberg, in S ources and Detecti on of Dark Matter in the Universe, edited by D.B. Cline (Springer Verlag, New York, 



to be published), astro- ph/0005265 



A.P Billyard and AA. Coley, Phys. Rev. D 61, 083503 (2000). 

L. Amendol a and D. Tocchini-Valentini, "Stationary dark energy: the present universe as a global attractor," astro- 



ph/0011243 



A. Burkert and J. Silk, in Dark Matter in Astro and Particle Physics, edited by H.V. Klapdor-Kleingrothaus and L. Baudis 
(IOP, Bristol, 1999). 

W.B. Lin, D.H. Huang, X. Zhang, and R. Brandenberger, Phys. Rev. Lett. 86, 954 (2001). 

D.N. Spergel and P.J. Steinhardt, Phys. Rev. L ett. 84, 3760 (200 0); J. P. Ostriker, ibid. 84, 5258 (2000); S. Hannestad, 



"Galactic halos of self-interacting dark matter," astro-ph/9912558 C. Firmani et al, Mon. Not. R. Astron. Soc. 315, L29 
(2000). 

J. Goodman, New Astronomy 5, 103 (2000). 

M. Kaplinghat, L. Knox, and M.S. Turner, Ph ys. Rev. Lett. 85, 3335 (2000). 



.vs. 

R. Cen, Astrophys. J. Lett, (to be published), p,stro-ph /0005206 



D. Pavon, and W. Zimdahl, Physics Letters A 179, 261 (1993). 
M.A.K. Gross, R.S. Somerville, J.R. Primack, J. Holtzman and A. Klypin, Mon. Not. R. Astr. Soc. 301, 81 (1998). 
W. Zimdahl, Mon. Not. R. Astr. Soc. 280, 1236 (1996). 

P. Blasi, in International Workshop on Observing Ultra High Energy Cosmic Rays from Space and on Earth, edited by 
(The American Institute of Physics, to be published). J. McDonald, "Non-Thermal WIMP Dark Ma tter, High Energy 



Cosmic Rays and Late-Decaying Particles From Inflationary Quantum Fluctuations," hep-ph/0011198 
D.J.H. Chung, E.W. Kolb, and A. Riotto, Phys. Rev. Lett. 81, 4048 (1998). 
V.S. Berezinsky, M. Kachelriess, and A. Vilenkin, Phys. Rev. Lett. 79 4302 (1997). 
VA. Kuzmin and VA. Rubakov, Phys. Atom. Nucl. 61, 1028 (1998). 
VA. Kuzmin and I.I. Tkachev, Phys. Rep. 320, 199 (1999). 
M. Birkel, and S. Sarkar, Astropart. Phys. 9, 297 (1998). 
V.S. Berezinsky, P. Blasi, and A. Vilenkin, Phys. Rev. D 58, 103515 (1998). 
P. Blasi, Phys. Rev. D 60, 023514 (1999). 



H. Ziaeepour, "Cosmic Equation of State, Quintessence and Decaying Dark Matter," astro-ph/000240C 



W. Zimdah l, D.J. Schwarz, A.B. Balakin, and D. Pavon, "Cosmic anti-friction and accelerated expansion," astro- 



ph/0009353 



W. Israel and W. Stewart, Ann. Phys. 118, 341 (1979); D. Pavon, D. Jou and J. Casas- Vazquez, Ann. Inst. H. Poincare 
A 36, 79 (1982); V. Romano and D. Pavon, Phys. Rev. D 50, 2572 (1994). 

R. Maartens, in Hanno Rund Workshop on Relativity and Thermodynamics, edited by S.D. Maharaj (Natal University, 
Durban, 1997). 



M. Anile, D. Pavon and V. Romano, "The case for hyperbolic theories of dissipation in relativistic fluids", gr-qc/9810014; 
D. Pavon, in Relativity and Gravitation in General, edited by J. Martin, E. Ruiz, F. Atrio, and A. Molina (World Scientific, 
Singapore, 1999). 

D. Pavon, D. Jou, and J. Casas- Vazquez, J. Phys. A 16, 775 (1983). 

WA. Hiscock and L. Lindblom, Ann. Phys. (N.Y.) 151, 466 (1983); Contemp. Math. 71, 181 (1988). 



12 



[38] G.F.R. Ellis and J. Wainwright in Dynamical Systems in Cosmology, edited by J. Wainwright and G.F.R. Ellis (Cambridge 
University Press, Cambridge, England, 1997). 

L.P. Chimento and A.S. Jakubi, Class. Quantum Grav. 10, 2047 (1993). 

L.P. Chimento and A.S. Jakubi, Class. Quantum Grav. 14, 1811 (1997); L.P. Chimento, A.S. Jakubi, and V. Mendez, Int. 
J. Mod. Phys. D 14, 1363 (1996); L.P. Chimento, A.S. Jakubi, V. Mendez, and R. Maartens, Class. Quantum Grav. 14, 
3363 (1997). 

A. Burd and A. Coley, Class. Quantum Grav. 11, 83 (1994). 
A. A. Coley and R.J. van den Hogen, Class. Quantum Grav. 12, 1977 (1995). 
S. Perlmutter et ai, Astrophys. J. 517, 565 (1999). 
A.G. Riess et ai, Astron. J. 116, 1009 (1998). 

G. Efstathiou, S.L. Bridle, A.N. Lasenby, M.P Hobson, and R.S. Ellis, Mon. Not. R. Astron. Soc. 303 L47 (1999). 
S. Perlmutter, M.S. Turner and M. White, Phys. Rev. Lett. 83, 670 (1999). 
L. Wang, R. Caldwell, J. P. Ostriker, and P.J. Steinhardt, Astrophys. J. 530, 17 (2000). 

P.J.E. Peebles, Principles of Physical Cosmology (Princeton University Press, Princeton, New Jersey, 1993). 
I. Maor, R. Brustein, and P.J. Steinhardt, Phys. Rev. Lett. 86, 6 (2001). 



P. Astier, "Can luminosity distance measurements probe the equation of state of d ark energy," astro-ph / 0008306 



M. Chevallier and D. Polarski, "Accelerated Universes With Scaling Dark Matter," |gr-qc/0009008 . 
J. A.S. Lima and J.S. Alcaniz, Astron. Astrophys. 348 1 (1999); M. Kaplinghat, G. Steigman, I. Tkachev, and T.P. Walker 
Phys. Rev. D 59 043514 (1999). 

M. Hamuy, M.M. Phillips, J. Maza, N.B. SuntzefT, R.A. Schommer, and R. Aviles, Astrophys J. 112, 2391 (1996). 

R. Lupton, Statistics in Theory and Practice (Princeton University Press, Princeton, 1993). 

P.M. Garnavich et ai, Astrophys. J. 509, 74 (1998). 

N. Kaiser, Astrophys. J. 498, 26 (1998). 

G.F.R. Ellis and M. Bruni, Phys. Rev. D 40, 1804 (1989). 

J. Bardeen, Phys. Rev. D 22, 1882 (1980). 

M. Bruni, P.K.S. Dunsby and G.F.R. Ellis , Astrophys. J. 3 95, 54 (1992). 
R. Maartens and J. Triginer, Phys. Rev. D 56, 4640 (1997); |gr-qc/9707018. 



13 



Am 



0.5 







-0.5 








0.2 



0.4 0.6 0.8 



FIG. 1. The magnitude residuals Am of the SNela (fit C data set of Ref. ]43|) from the best fit superattractor QDDM model 
with {/3, M) = (0.395, 23.96). For comparison is shown the magnitude residual for the best fit ACDM model, corresponding to 
(n m o,^A,M) = (0.54,1.09,23.93) @. 
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FIG. 2. The estimated probability density distribution (normalized likelihood) for the acceleration parameter f3 of 
superattractor model. 
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